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We report the first global study of dynamic correlation energies (DCEs) associated with rotational 
motion and quadrupole shape vibrational motion in a covariant energy density functional (CEDF) 
for 575 even-even nuclei with proton numbers ranging from Z = 8 to Z = 108 by solving a five¬ 
dimensional collective Hamiltonian, the collective parameters of which are determined from triaxial 
relativistic mean-field plus BCS calculation using the PC-PKl force. After taking into account these 
beyond mean-field DCEs, the root-mean-square (rms) deviation with respect to nuclear masses is 
reduced significantly down to 1.14 MeV, which is smaller than those of other successful CEDEs: 
NL3* (2.96 MeV), DD-ME2 (2.39 MeV), DD-ME(5 (2.29 MeV) and DD-PCl (2.01 MeV). Moreover, 
the rms deviation for two-nucleon separation energies is reduced by ~ 34% in comparison with 
cranking prescription. 
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Nuclear energy density functional theory is nowadays 
one of the most important microscopic approaches for 
large-scale nuclear structure calculations in medium and 
heavy nuclei based on a universal energy density func¬ 
tional (EDF) with a few parameters constrained by the 
properties of finite nuclei and nuclear matter or neutron 
stars. It stands out as a unique microscopic model that 
can describe not only the masses of all existing nuclei, 
but also other key quantities for simulating nucleonthe- 
sis process, including beta decay rate and fission rate in 
a unified way [1]. However, it is still a challenge for the 
current implementation of EDF to achieve satisfied accu¬ 
racy. Therefore, great efforts have been devoted in many 
aspects to improving the accuracy of the EDF for atomic 
nuclei and subsequently deepening our understanding on 
the origin of elements in universe. 

Nuclear binding energy or mass is one of the most fun¬ 
damental properties of atomic nuclei. The root-mean- 
square (rms) deviation with respect to the measured nu¬ 
clear masses in the EDF is typically around several MeV. 
Only after taking into account the beyond mean-field dy¬ 
namic correlation energies (DCEs) in global fitting of the 
EDF to nuclear masses can one achieve the rms deviation 
of a few hundred keV [2, 3] by keeping a good descrip¬ 
tion of nuclear matter properties. In these studies [2, 3], 
the DCEs related to rotational and vibrational motions 
have been included phenomenologically with the crank¬ 
ing prescription. A better treatment of these DCEs is to 
carry out calculation with exact quantum number projec¬ 
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tion and generator coordinate method (GCM), which has 
been done based on a Skyrme SLy4 force [4, 5] with the 
assumption of topological Gaussian overlap approxima¬ 
tion (GOA) or based on Gogny DIM and DIS forces [6] 
without using the GOA. In both the beyond mean-field 
(BMF) calculations, only axially deformed configurations 
are included and the obtained quadrupole DGEs are on 
average 3 — 4 MeV. Although the BMF approaches with 
exact projections and GGM for triaxially deformed nu¬ 
clei have already been developed in recent years [7-9], 
they cannot be adopted for a large-scale calculation of 
the DGEs because of extreme time-consuming. In con¬ 
trast, the method of five-dimensional collective Hamilto¬ 
nian (5DGH) with parameters determined by the mean- 
field calculations is much cheaper in numerical realization 
and this method has been adopted to study low-lying 
states for a large set of even-even nuclei based on the 
Gogny DIS force [10] and to optimize the new Gogny 
DIM force with a great success [11]. 

Besides the nonrelativistic EDF, the covariant energy 
density functional (GEDF) has received wide attention 
and achieved great success in describing many phenom¬ 
ena of both stable and exotic nuclei [12-15]. However, 
the description of nuclear binding energy based on the 
popular GEDFs is still not satisfactory. It has been 
found in the most recent large-scale calculation of nuclear 
binding energy with axially deformed relativistic Hartree- 
Bogoliubov (RHB) for all proton number Z < 104 even- 
even nuclei that the rms deviations with respect to the 
640 measured masses of even-even nuclei in the AME2012 
mass evaluation are 2.96 MeV for NL3*, 2.39 MeV for 
DD-ME2 and 2.29 MeV for DD-ME^, and 2.01 MeV for 
DD-PGl [16]. Similar accuracy has been found in other 
global studies of nuclear masses with axially deformed 
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RMF+BCS approaches [17-19]. 

Therefore, a natural question to be asked is: to which 
extent the DCEs can improve the description of nuclear 
masses based on the existing CEDFs. To address this 
question, we have recently carried out a global calcula¬ 
tion of the DCEs for the 575 even-even nuclei with axially 
deformed relativistic mean-field plus BCS (RMF-^BCS) 
approach using the point-coupling PC-PKl parametriza- 
tion [20]. The DCEs of deformed nuclei have been eval¬ 
uated at the mean-field level with the cranking prescrip¬ 
tion adopted in building the mass tables based on Skyrme 
EDFs [ 2 , 3]. After including these correlation energies, 
the rms deviation of the masses is reduced from 2.58 MeV 
to 1.24 MeV [20]. However, it should be pointed out that 
the cranking prescription for the DCE based on mean- 
field ground-state wave function is improper for the nu¬ 
clei with transitional characters and for the nuclei with 
other energy minima competing with the ground-state 
one. Because the resultant DCE correction to binding 
energy in these cases depends much on the choice of the 
ground-state configuration which is, however, ill-defined 
and could be easily altered by the DCE [20, 21 ]. More¬ 
over, for the transitional nuclei, large shape-fluctuation 
effect on nuclear properties is expected. All these defi¬ 
ciencies in the cranking prescription can in principle be 
taken care by the 5DCH method which is regarded as a 
GOA of angular momentum projection plus GCM. The 
5DCH with the parameters determined by the CEDE has 
been implemented in Ref. [ 22 ] and the success of it has 
been illustrated in a series of calculations for the spheri¬ 
cal, transitional, and deformed nuclei from A ^ 40 to su¬ 
perheavy region [ 21 , 23-28]. Recently, we have performed 
a detailed comparison between the 5DCH and exact pro¬ 
jections plus GCM calculations for the low-lying states 
of ^^Kr based on the PC-PKl force and found that these 
two methods gave similar results [29]. 

In this work, we revisit the DCEs in the CEDE us¬ 
ing the 5DCH method for the 575 even-even nuclei with 
8 < Z < 108 and examine the influence of these BMF cor¬ 
relation energies on the predicted masses and separation 
energies. Compared with our previous work on the DCEs 
based on the cranking prescription [ 20 ] , the present work 
takes into account more correlations in a more proper way 
as mentioned above and this paper presents the most ad¬ 
vanced systematic calculation of the DCEs in the CEDE. 

To this end, we first carry out a large-scale deformation 
constrained triaxial RMF+BCS calculation to generate 
the mean-field wave functions in the whole (/3, 7 ) plane, 
where the PC-PKl force [30] in the particle-hole chan¬ 
nel and a separable paring force [31-34] in the particle- 
particle channel are adopted. This separable pairing 
force has finite range and preserves translational invari¬ 
ance. In particular, the numerical calculation is much 
simpler than that of Gogny force, which makes it more 
suitable for systematic studies. Several studies have al¬ 
ready shown that the separable pairing force is quite suc¬ 
cessful in the description of nuclear shape transition and 
low-lying spectrum [ 21 , 25, 26]. The Dirac equation for 


single-particle wave function is solved by expanding in a 
set of eigenfunctions of 3D harmonic oscillator in carte¬ 
sian coordinate with 12, 14, and 16 major shells for nu¬ 
clei with Z < 20 , 20 < Z < 82, and Z > 82, respectively. 
The single-particle wave functions, occupation probabil¬ 
ities, and quasiparticle energies are used to calculate the 
mass parameters, moments of inertia, and collective po¬ 
tentials in the 5DCH, all of which are functions of the 
deformations /3 and 7 . The properties of nuclear ground 
state Oj^ and other low-lying excited states are obtained 
by solving the 5DCH [ 22 , 23]. The DCE Ecorr is defined 
as the energy difference between the mean-field ground 
state and the 5DCH ground state Consid¬ 

ering the facts that the PC-PKl force was parameterized 
to many spherical nuclei ranging from O to Pb and the 
5DCH method for nuclei around double-closed shell is 
improper, these lead to erroneous negative DCEs. For 
these nuclei, the Ecorr is imposed to be zero in this work. 
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FIG. 1: (Color online) Contour map of the quadrupole dy¬ 
namical correlation energies Ecorr calculated by the CEDF 
based 5DCH model as functions of the neutron and proton 
numbers. 


Figure 1 displays the contour map of the quadrupole 
dynamical correlation energies Ecorr calculated by 5DCH 
based on CEDF with PC-PKl force. The quadrupole dy¬ 
namical correlation energies range from 0 to ^ 4.4 MeV, 
and mostly vary in the region of 2.0 ^ 3.5 MeV. For 
the semi-magic nuclei with Z = 28, 50, 82 and V = 28, 
the correlation energies are nonzero or even rather large. 
This is because the energy surfaces of these nuclei are 
either soft around energy minimum or with shape coex¬ 
isting phenomena. The Ecorr of the transitional nuclei 
with Z ^ 54, 78 and N ^ 60, 90 is pronounced, generally 
larger than 2.5 MeV due to the shape fluctuation. The 
Ecorr of the well-deformed nuclei with A ^ 170 is around 
2.5 MeV and it reduces to ~ 2.0 MeV for the heavier 
nuclei with A r\j 250. 

In Fig. 2, we compare our dynamical correlation en¬ 
ergies with those calculated from 5DCH model based on 
Hartree-Fock-Bogoliubov (HFB) calculation using Gogny 
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FIG. 2: (Color online) The quadrupole dynamical correlation 
energies calculated by 5DCH based on CEDF with PC-PKl 
force (open circles) in comparison with those based on HFB 
using Gogny DIS force (upper triangles). 


DIS force [10]. The systematics of E'corr are similar in 
both calculations. However, the DCEs from PC-PKl 
force are systematically smaller than those from Gogny 
DIS force, and the rms deviation of these two results is 
^ 1.95 MeV. The difference in the DCE might be origi¬ 
nated from the collective parameters (especially the zero 
point energies), which are sensitive to the effective inter¬ 
actions, in particular, to the pairing properties [35]. 
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EIG. 3: (Color online) Discrepancy of the CEDE calculated 
binding energies by PC-PKl with the data for 575 even-even 
nuclei [36]. In panel (a), the CEDE calculated binding ener¬ 
gies are given by the binding energies of the lowest mean-field 
states, while in panel (b) the dynamical correlation energies 
are taken into account. 


Figure 3 displays the discrepancy of the calculated 
binding energies with respect to the data for 575 even- 
even nuclei [36] . The mean-field binding energy (denoted 
as Hxri) is obtained by adding the energy correction from 
static triaxial deformation to that from our previous axi¬ 
ally deformed RMF+BCS calculations [20]. Such energy 
correction is defined as the energy difference between the 
ground state and the lowest non-triaxial minimum of the 
triaxial deformed nucleus (totally 31 in the present work) 
calculated by the triaxial RMF+BCS with the same in¬ 
puts as in Ref. [20] . The static triaxial deformation effect 
makes the nucleus more binding in an average value of 
^0.36 MeV, closer to the experimental data. Our main 
findings are as follows: 

(i) The mean-field binding energies are systematically 
underestimated by about 3 MeV. By including the 
DCEs from the 5DCH calculation as shown in 
Fig. 1, the rms deviation is reduced from 2.52 MeV 
to 1.14 MeV, which is about 100 keV improvement 
in comparison with the cranking prescription in 
Ref. [20]. The remaining discrepancy between the 
calculated and measured binding energies for most 
nuclei is in between —1.0 and 1.0 MeV. 

(ii) The large discrepancy for nuclei in the neighbor¬ 
hood of (Z, N) ^ (40,40) in the cranking approx¬ 
imation calculation [20] has been modified signifi¬ 
cantly in the present calculation. This is because 
the shape coexistence and large shape fluctuation in 
this mass region are considered in 5DCH method. 

(iii) Similar as in Skyrme-HFB [2, 3] and Gogny- 
HFB [11] mass models, large deviations occur 
around magic numbers. In particular, the bind¬ 
ing energies of V ^ 126 isotones are significantly 
overestimated, which has also been found in other 
CEDFs [16]. It is still an open question which has 
not been understood. 

Figure 4 shows the discrepancy of the CEDF calcu¬ 
lated two-neutron and two-proton separation energies by 
PC-PKl with respect to the data [36]. It is remarkable 
that the separation energies are reproduced very well ex¬ 
cept for the nuclei around N ^ 20,126 and Z ^ 28,82. 
The discrepancy between the calculated and measured 
values for most nuclei is in between -0.5 and 0.5 MeV, 
and the rms deviations are 0.62 MeV and 0.69 MeV for 
S 2 n and 5'2p, respectively, which are smaller than the 
values 0.96 MeV and 1.03 MeV in Ref. [20] based on 
the cranking prescription. For comparison, we also com¬ 
pute the separation energies from the binding energies of 
the lowest mean-field states (i.e. B^ri in Fig. 3) and the 
rms deviations are respectively 0.91 MeV and 0.96 MeV 
for neutrons and protons. In short, the DCEs from the 
5DCH calculation gives a significant improvement for the 
description of two-nucleon separation energies. 

In summary, we have carried out a global study of 
the dynamical correlation energies associated with the 
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FIG. 4: (Color online) Discrepancy of the theoretical two- 
neutron (panel a) and two-proton (panel b) separation ener¬ 
gies extracted from the calculated binding energies including 
the DCEs with respect to the data [36]. 


quadrupole shapes for the 575 even-even nuclei with pro¬ 
ton numbers ranging from Z = 8 to Z = 108 by solv¬ 
ing the 5DCH with the collective parameters determined 
from the CEDF calculation using PC-PKl force. The 
DCEs range from 0 to ^ 4.4 MeV, and mostly vary in 
the region of 2.0 ^ 3.5 MeV. After including DCEs, the 
CEDF predictions for 575 masses, 521 two-neutron sep¬ 
aration energies, and 497 two-proton separation energies 
are improved significantly with the rms deviations 1.14 
MeV, 0.62 MeV, and 0.69 MeV, respectively. It is re¬ 
markable that the rms deviation for two-nucleon separa¬ 
tion energies is reduced by ^ 34% in comparison with 
our previous results based on cranking prescription. 
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